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MATLAB 简 介 

信号 的 MATLAB 表 示 

利用 MATLAB 进 行 系统 的 时 域 分 析 
利用 MATLAB 进 行 信号 的 频 域 分 析 
利用 MATLAB 分 析 系 统 的 频率 特性 
利用 MATLAB 进 行 连续 系统 的 s 域 分 析 




















利用 MATLAB 进 行 离散 系统 的 z 域 分 析 
利用 MATLAB 进 行 系统 的 状态 变量 分 析 
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MATLAB 简 介 
(Matrix Laboratory ) 
MATLAB 的 工作 方式 
如 何 获取 帮助 
表达 式 一 一 变量 、 数 值 、 算 数 运算 符 、 
关系 运算 符 、 示 辑 运 算 符 、 冒 号 运算 符 
数组 及 其 运算 
函数 文件 
循环 (FOR、WHILE 循环 ) 
基本 绘图 语 铝 





一 、MATLAB 的 工作 方式 








cf) 窗 日 命令 方式 
《2) 运行 以 .M 为 扩展 名 磁盘 文件 





下 作 方 式 。 罗 


多用 plot 函 数 画 一 个 方 波 z 
t=[-100113]; 
Xx=[001100]; ni 1 
plot(bX); 

Xlabel(t);ylabel(x(b); 

axis([-1302]); 





二 、 获 取 帮 助 











命令 窗口 输入 : ”help+ 函 数 名 
h 


例如 help plot 











养 变量 

不 需要 变量 的 类 型 说 明 

变量 名 的 第 一 个 字符 必须 是 字母 

w 变量 名 长 度 : 不 超过 31 个 字符 

w 大 写 和 小 写 的 字母 视 为 不 同 的 字符 
例如 : num_students = 25 

特殊 变量 : 
pi 表示 圆周 率 ，inf 表示 无 穷 大 ，NaN(Not a 
NumbenD 表 示 不 定量 ， 如 0/0。 

















数值 全 
MATILAB 用 常规 的 十 进 制 表示 数值 


”用 i 或 j 作 为 后 缀 来 表示 复数 的 虚 部 
例 1.23$e5 表 示 1.235Sx105，Xx=2+3j 


abSs(X) 求 复数 x 的 模 

angle(X) 求 复数 x 的 相 角 ( 驳 度 ) 
real(X) 求 复数 x 的 实 部 
Imag(X) 求 复数 x 的 虚 部 
conj(x) 求 复数 x 的 共 斩 





三 、 表 达 式 








姜 运算 符号 
> 算数 运算 符 

依 手 加 

4 - 城 

多 亚 乘 

4 / 除 

4 ^ 乘 方 

4 “′ 和 拢 阵 的 复 共 斩 转 置 








三 、 表 达 式 


养 运算 符号 
> 逻辑 运算 符 
4 AK&B 逮 辑 与 and) 
4 AIB 逻辑 或 (on 
4 ~A 逻辑 非 oob 











三 、 表 达 式 








状 运 算 符号 
> 关系 运算 符 

多 A<B 小 于 

4 二 人 
ES 
> 了 三 二 并 兰 于 
9 A== 等 于 

4 A-=B 不 等 于 











养 运算 符号 
> 冒号 运算 符 
4 表达 式 1:10 表示 产生 一 个 行 向 量 ， 它 的 值 为 
Eee 
4 表达 式 10:-2:1 表示 产生 一 个 递减 的 行 向 量 ， 
它 的 值 为 
10 8 6 4 2 











四 、 数 组 








芋 ， 数 组 的 构造 


例 =x=2:5= 产生 一 个 数组 ， 它 的 值 为 
X(D)=2,X(2)=3;X(3)=4;X(4)=5 


X=linspace(0,2,11) 将 区 间 [0，J] 均 匀 抽 样 11 点 作为 数组 x 
给 2 维 数组 赋值 时 ， 用 分 号 表示 一 行 的 结束 ， 
zEEEZE 人 E23E 人 和 





四 、 数 组 








1.， 数组 的 构造 
MATLAB 提供 了 王 些 产生 基本 和 矩阵 的 函数 
ZerOS 产生 年 阵 元 素 全 为 0 的 玫 阵 
ones 产生 珑 阵 元 素 全 为 1 的 和 矩 阵 
rand 产生 (0,1) 均 匀 分 布 随机 数 拖 阵 
randn ”产生 正 态 分 布 随机 数 和 矩阵 











四 、 数 组 








2. 数组 的 运算 
数组 和 一 个 标量 相 加 或 相 乘 
例  y=X-1 Z=3#X 
2 个 数组 的 对 应 元 素 相 乘除 . 涯 ./ 
例 “ z=X.*y 
确定 数组 大 小 的 函数 
Size(A) 返回 值 数组 A 的 行 数 和 列 数 三维) 
length(B) 确定 数组 了 B 的 元 素 个 数 一 维 ) 





五 、 函 数 文件 








二 M 文 件 的 第 一 行 包含 function 
功能 : 建立 一 个 函数 ， 可 以 同 MATLAB 的 库 
函数 一 样 使 用 。 





五 、 函 数 文件 








例 ， 编 一 个 绘制 图 示 波 形 的 函数 。 小 
1 0 1 


调用 函 数 ti， 并 画 出 它 的 波形 





六 、For 循环 








例 : 编写 计算 s=1+2+3+...+100 的 MATLAB 程 序 


S=0; 





七 、While 循环 








例 : 计算 *= 瑟 十 吉 + 到 +… 的 值 ， 目 误差 小 于 104 
S=f 人 
D=1]; 
eps=1le--0; 





fprintt(Cs=%.Sfn',S) 








八 、plot 函 数 绘图 函数 (continuous) 








t=]linspace(0,4*pl312); 
plot(tsin(b ,bcos(D，-.); 
title(my figure ); 
Xlabel(Ct); 
legend(Csin(b' ,cos(D); 





绘图 函数 (continuous) 





八 、plot 函 数 






































my fure 
TS 
Sin 人 
= = 一 《0s 人 
! - 
1 
7 8 ! 
7 
7 
-0.z 1 《 7 
-0.4 1 1 
-0E 1 A f 
-0E A 1 
1 RN 
0 区 8 加 12 14 


九 、stem 函 数 





绘图 函数 (discrete) 








k=0:39; 
Stem(k,cos(0.9*#pixk)); 
title(Ccos(0.9\pik)); 





可 
所 





九 














cos(0.9zk) 


0.8 上 
0.6 上 
0.4 上 
0.2 上 

0 上 


-0.2 上 1 | 
-0.4 上 
-0.6 上 














站 信号 的 MATLAB 表 示 


4 基本 信和 号 的 MATLAB 表 示 
指数 信号 4ew 、 指 数 序列 履 、 抽 样 函数 Sa(D、 
正 弱 型 信号 、 拢 形 脉冲 信号 、 三 角 脉 冲 信号 


4 信号 基本 运算 的 MATLAB 实 现 
尺度 变换 、 翻 转 、 时 移 、 相 加 、 相 乘 、 
差分 与 求 和 、 微 分 与 积分 





一 、 基 本 信号 的 MATLAB 表 示 





> 指数 信号 Aew y= 三 Atkexp(akft); 
> 指数 序列 圭 运 算 a.Ak 实 现 
> 正弦 型 信号 内 部 函数 cos() 和 sin() 


> 抽样 函数 Sa(D) sinc( 

> 天 形 脉冲 信号 y=ftectpuls(twidthb) 
> 三 角 波 脉冲 信和 号 

y=trpuls(t width,skew) 








一 、 基 本 信号 的 MATLAB 表 示 





%decaying exponential 
t=0:001:10; 
A=1; 

a=--0.4; 
ft=A*xexp(axt); 
plot(tfb 
t=0:0.1:10; 
A=1; 

a=--0.4; 
ft=Ax*xexp(axt); 
stem(t,fb) 











[和 

















[ 关 汪 Te 


一 、 基 本 信号 的 MATLAB 表 示 











% rectpuls 
t=0:0.001:4; 





械 三 此 
ft=rectpuls(t 一 2#TT); 














plot(tfbD 








axis([0,4,-0.5,1.5]) 


一 、 基 本 信号 的 MATLAB 表 示 














2% tripuls 司 

人 -3:0.001:3; o 
ft=tripuls(t,4.0.5); 
plot(t,fb) | 到 





ft=tripuls(t,4,1); 











基本 信号 的 MATLAB 表 示 











色 unit impuls sequence 
k=-50:530; 
delta=[zeros(1,90),1,zeros(1,50)]; 
stem(k,delta) 





function [人 kj=impseq(kOKkl,k2) 
% 产 生 frklj=delta(k-k0); 站 
k=[Kkl:k2]; 仁 [(k-kO)==0]; 
k0=0:k1=-S0;:k2=50; 
[=impseq(kO,.K1,k2); 

stem(k, 全 


一 、 基 本 信号 的 MATLAB 表 示 








% unit step Sequence 
k=-50:50; 

UKk=[zeros(1,390), ones(1,31)]; 
stem(k:uk) 





function [人 kj]=stepseq(kOkl,k2) 
包产 生 ffkl]=u(k-k0);kl1<=k<=k2 
k=[kl:k2]:f=[(k-kO0)>=0]; 
k0=0:k1=-S0:K2=50; 
[fk]=stepseq(kOkl,k2); 
stem(k, 旭 


二 、 信 号 基本 运算 的 MATLAB 实 现 








1， 信和 号 的 尺度 变换 、 翻 转 、 时 移 〈 平 移 ) 
全 -3:0.001:3; 
从 1=tripuls(2#t4.0.5); 
Subplot(2,1,1) 
plot(tftl) 
title(Xx(207) 
ft2=tripuls((2 一 2*#0,4.0.5); 
Subplot(2,1,2) 
plot(tft2) 
title( X(2-20) 


已 知 三 角 波 x(D， 用 MATLAB 画 出 的 xC2D 和 x(2-27 波形 


x(2h 








mohooecoo 


oa 品 
N ， 


Xx(2-2t) 




















borooocoo 


2 晤 
口 
口 
mm 
四 


二 、 信 号 基本 运算 的 MATLAB 实 现 








2.， 信 号 的 相 加 与 相 乘 
w 相 加 用 算术 运算 符 “+2” 实 现 
相 乘 用 数组 运算 符 “. 关 "实现 


例 : 画 信号 Aewcos(wW 夺 幼 的 波形 
t0:0.001:8; 
A=1l; a=--0.4; 
w0O=2*piphi=0; 
ft1=Ax*xexp(axt).*sin(WOx*t+phi); 
plot(tftl) 

















-0 
(了 


二 、 信 号 基本 运算 的 MATLAB 实 现 








3， 离 散 序列 的 差分 与 求 和 


连续 信号 的 微分 与 积分 
> 差分 y=diffd); 
羡 求 和 y=Sum(f(k1l:k2)); 


> 微分  y=diff(pyjh; hb 为 数值 计算 所 取 时 间 间 隔 
三 定 积分 quad('function_name”,a.b); 
function_name 为 被 积 函 数 名 ，a 和 b 指 定 积 分 区 间 。 





二 、 信 号 基本 运算 的 MATLAB 实 现 








3， 离 散 序列 的 差分 与 求 和 
连续 信号 的 微分 与 积分 


例 :已 知 三 角 波 xD， 画 出 其 微分 与 积分 的 波形 








%differentiation %integration 

h=0.001;t -3:h:3; 个 二 3.0L3; 

yl1=diff(f2_2(t))*1/h; for x=1l:length(b) 

plotkt(l:length(D-1D,yl) y2(x)=quad(f2_2, -3,t(Cx)); 
end 


plot(ty2) 


三 角 波 刀 D) 微 分 与 积分 的 波形 























0 dx(bDydt 

5 r + r + + 

0 三 = 国 
0.5 

j| 
15 。 

四 7 1 0 公 3 
入 integral of x(U 
5 - 

















生 利用 MATLAB 进 行 系统 的 时 域 分 析 


* 连续 时 间 系 统 零 状态 响应 的 求解 

4 连续 时 间 系 统 冲 激 响应 和 阶 跃 响应 的 求解 
4 离散 时 间 系 统 零 状 态 响应 的 求解 

* 离散 时 间 系 统 单位 脉冲 响应 的 求解 

4 离散 卷 积 的 计算 


一 、 连 续 时 间 系 统 雯 状态 响应 的 求解 





lsim 





t 表示 计算 系统 响应 的 抽样 点 向 量 
x 是 系统 输入 信号 向 量 ， 
sys 是 LTI 系 统 模型 ， 借 助 tt 函数 获得 





b 和 a 分 别 为 微分 方程 右 端 和 左 端 各 项 的 系数 向 量 


03y (+Tay (DT+TwWy (GD)+aoyG) 三 DT (DPCD+TDX (CD+DPox(GD 





二 、 连 续 系统 冲 激 响应 和 阶 跃 响应 求解 








连续 时 间 系 统 冲 激 响应 可 用 impulse 函 数 直接 求 出 ， 
其 调用 形式 为 


Impujse 





连续 时 间 系 统 阶 跃 响应 可 用 step 函 数 直接 求 出 ， 
其 调用 形式 为 





t 表示 计算 系统 响应 的 抽样 点 向 量 
SYS 是 LTI 系 统 模型 





三 、 离 散 时间 系 统 雯 状态 响应 的 求解 








b ,a 分 别 是 差分 方程 左 、 右 端的 系数 和 问 量 
x 表示 输入 序列 ，y 表示 输出 序列 
羡 woK-= 立 5 


可 用 MATLAB 表 示 为 





四 、 离 散 时 间 系 统 单 位 脉冲 响应 的 求解 








b, a 分 别 是 差分 方程 左 、 右 端的 系数 向 量 
k 表示 输出 序列 的 取 值 范围 
h 就 是 单位 脉冲 啊 应 





五 、 离 散 卷 积 的 计算 


式 中 ab 为 竺 卷 积 两 序列 的 向 量 表示 ，e 是 卷 积 结果 。 


例 : (S3+2s+3)(S2+3s+2) 可 用 下 面 MATLAB 语 名 求 出 
a =[1.0.2.3]; 
b =[1.3.2]; 
c=conyv(a,b) 











例 1 求 系统 >"(D+2y(D+100y(D=10x(D 的 零 状 态 
响应 ， 已 知 x(D=sin(2TD 2(D。 


%program3_1 微 分 方程 求解 
ts=0:te=95:dt=0.01; 光 
SySs=tt([10],[1 2 100]); 
t=ts:dt:te; 
X=Sin(2*#plxt); 
y=lSim(SySs,X,D; 
Plot(by); 
Xlabel(CTime(sec)) 
ylabel(y(O) 








2 3 
Timelsec) 


例 2 求 系统 y" (D)+2y' (D+100y(O=IOx(D) 的 零 状 
态 啊 应 ， 已 知 x(D =6GD 。 


%program3_2 连 续 时 间 系 统 的 冲 激 响应 
ts=0:te=9:dt=0.01; 
SySs=tft([10],[1 2 100]); 








t=ts:dt:te; 史 
y=lmpulse(SySs,D); Hi 
plot(by); 

xlabel(Time(sec)) 业 
ylabel(h(b) 


Timelsec) 


例 3 分 析 噪 声 干 扰 的 信号 x[ 间 奈 [ 同 +d[ 各 通过 M 点 滑动 
平均 系统 的 响应 ， 7 = 号 二 2 人 
其 中 s[ 癌 =(290.9* 是 原始 信号 ，d[ 问 是 噪声 。 


R =31 ;d=rand(],R) -0.3; 
k=0:R 一 1; 

S=2#K.#(0.9.^k); X=S+d; 

皇 gure(]); 
plot(k,d,T-. ,Sb-- 上 ,X, 8-); 
M =5;b = ones(M,1)/M; a= 1; 
y=filter(b,a,x); 

全 gure(2); 

plot(k,s,b-- ,kk,y'T-); 


例 3 分 析 噪 声 干 扰 的 信和 号 8 量 过 M 点 滑动 
平均 系统 的 响应 ， xD= 寺 忆 x[ 大 一 吕 ] 
其 中 s[ 辣 =(20.9* 是 原 始 信号 ” d[ 问 是 噪声 。 


噪声 干扰 信号 x[ 如 =s[ 轨 +d[ 回 通过 M 点 滑动 平均 系统 的 响应 















































20 30 20 30 
Timeindexk Timeindex k 


例 4 求 系统 y[K]+3y[ 人 1]+27[1]=10x[ 如 的 单位 











肪 证 啊 应 。 
% program 3_4 离散 系统 的 单位 脉冲 响应 
k=0:10; 单位 及? 
a=-[1 3 2]; 3000 天 一 
b=[10]; 2000 
Ph=impz(b,ak); 1000 
stem(k,b) ULUREOOOO 
-100 
-200 


例 $S 计算 x[A* y[ 则 并 画 出 卷 积 结果 ， 已 知 
X[N={1.2.3,4; 寺 0.12.3}， 
y[6={HL1111; 季 0,1.2.3.4] 。 





% program 3_5 
X=[1;2,3,4]; 
YE[ULSLSLSISU 
Z=CcCOnV(X,y); 
N=length(Z); 

















stem(O:N 一 1,Z); 


半 利用 MATLAB 进 行 信号 的 频 域 分 析 


* 周期 信号 频谱 的 MATLAB 实 现 
* 用 数值 积分 分 析 非 周期 信号 频谱 


一 、 周 期 信号 频谱 的 MATLAB 实 现 








频谱 Fn 一 般 为 复数 ， 可 分 别 利用 abs 和 angle 函 
数 获得 其 幅 频 特 性 和 相 频 特性 。 


其 调用 格式 分 别 为 


周期 信号 的 频谱 Cn 为 离散 信号 ， 可 以 用 stem 
8 其 频谱 图 。 





一 





例 1 试用 MATLAB 画 








审 号 的 频谱 。 








画 三 角 波 信号 频谱 的 MATLAB 程 序 


N=8; 

nl= -N:-1; 和 计算 n=-N 到 -1 的 Fourier 系 数 
cl= 一 43j*sin(n1*#pi/2)/piA2./n1.^2; 

c0=0; 匈 计 算 n=0 时 的 Fourier 系 数 

n2=1:N; % 计 算 n=1 到 N 的 Fourier 系 数 

c2= -4#jssinn2#pi/2)/piA2./n2.A2; 

cn=[cl c0 c2]; 

n= -N:N; 

Subplot(2,1,1); 

stem(nabs(cn));ylabel(CCn 的 幅度 ); 
Subplot(2,1,2); 

stem(n,angle(cn)); 

ylabel(CCn 的 相位 );xlabel(\omegaNAomega07); 

















0.5 


0.4 下 








0 


@/o0 


E 形 脉冲 的 Fourier 级 数 表示 式 。 并 用 
的 信 








例 2“ 求 周期 匀 
MATLAB 求 出 由 前 NM 项 Fourier 级 数 系数 得 昌 


号 近似 波形 。 









多 
2 


你 
C, = 二 Sa(twuz/12) 


(三 2 Ge 有 十 泛 sa 和 COSTLC 


1=] 





肥 帮 三 | 辐 7 三 2 天/ 三 区 
立 OD=0.5+》， SaC) cos(nrn) 


7m=1 


2 Gibbs phenomenon 


t=-2:0.001:2; 

N=input(Number of harmonics= ); 

cO=0.3; 

XN=c0*ones(1l,length(b)); %dc component 

forn=0:2:N ”% even harmonics are zero 
XN=XN+cos(pixn#xt)*#sinc(n/2); 

end 

plot(xN); 


2 Gibbs phenomenon 





中 


N=5 o N=15 








NE50 " 中 jV=500 



































二 、 用 数值 积分 分 析 非 周期 信号 频谱 








数值 函数 积分 quad8 可 用 来 计算 非 周期 信号 频谱 
。 y= quad8(F'ab) 


E 是 一 个 字符 串 ， 它 表示 被 积 函 数 的 文件 名 。 


ab 分 别 表示 定 积 分 的 下 限 和 上 限 
quad8 的 返回 是 用 自 适应 Simpson 算 法 得 出 的 积分 值 











例 3 ”试用 数值 方法 近似 计算 三 角 波 信号 的 频谱 





: 图 示 三 角 波 可 表示 为 
(=G= 负 挛 (D 
三 角 波 信和 号 频谱 的 理论 值 为 
XUJ= Sa2(w/12) 








例 3 ”试用 数值 方法 近似 计算 三 角 波 信号 的 频谱 


function y=sfl(tw); 
y=(t>=-1 久 t<=]1).*#(1-abs(t)).#xeXxp(-j*wWxt); 


w=linspace(-6#pi,6*pi,512); 
N=length(w);X=zeros(1,N); 
for k=1:N 

X(oO=quad8(Csfl ,-1,1, 口 , 吕 ,wdo)); 
end 
figure(l); 
plot(wreal(X));title() 
Xlabel(\omega);ylabel(XGU\omega)); 
人 gure(2); 
plot(wreal(X)-sinc(w/2/pi).^2); 
xlabel(\omega):title(" 计 算 误 差 ); 


X(jey 








例 3 ”试用 数值 方法 近似 计算 三 角 波 信号 的 频谱 


二 
运行 结果 

业 

8 
1 寺 5 5 5 志 瑟 及 
0 计算 误差 

ee 

5 

妹 

-20 15 10 吐 0 5 10 15 20 


站 利用 MATLAB 进 行 系统 频 域 分 析 


* 连续 系统 频率 啊 应 的 计算 
4 周期 信号 通过 系统 的 响应 
* 离散 系统 频率 响应 的 计算 


一 、 连 续 系统 频率 响应 的 计算 








B(O _ DDGw ”+DC2)(0j0O) ”+ …DCM +D 
4(O  a((jo)> +aC2)(Go > +…aCN+D 


计算 频 响 的 MATLAB 函 数 


互 UjOw) = 








b 分子 多 项 式 系数 

a 分 母 多 项 式 系数 

w 需 计 算 的 瑟 G 四 的 抽样 点 
(数组 w 中 少 需 包 含 两 个 ow 的 抽样 点 )。 





一 、 连 续 系统 频 响 特性 的 计算 








例 芋 三 阶 归 一 化 的 Butterworth 低 通 滤 波 器 的 系统 函数 为 
| 
互 UJ = 一 一 = 
WUO) +2UO) +20jO+1 

试 画 出 IEGwl 和 mg)。 
w=linspace(0,3,200); 
= 站 由 3 三 [加 2 
h=freqs(b,a,w); 
Subplot(2,1,1); 
plot(w,abs(h)); 
Subplot(2,1,2); 
plot(w,angle(h)); 





一 、 连 续 系统 频 响 特性 的 计算 





子 阶 Butterworth 低 通 滤波 器 的 幅度 响应 和 相位 响应 





IHGo)l 








9(o) 





@ (rad/S) 


二 、 周 期 信号 通过 系统 的 响应 
例 周期 方 波 通过 RC 系 统 的 响应 。 






































X(b) 人 

1 下 

qd) CE 
EEC IT 3 5 

态 河 LjaC 1 C， 2 1C0T sa 玫 

R+I/joc 1+jwRC 7 5) 有 5 
yD= 二 CEGnm)e" 05+SSa( 加 | -|cosozaaz+ 内 

I 王 - = 2 |1+jzRR 


二 、 周 期 信号 通过 系统 的 响应 








例 2 周期 方 波 通过 RC 系 统 的 响应 。 
%p5_2 Periodic Signal pass LII system 
T=4;w0=2*pi/T:RC=0.1; 
t= -0:0.01:6;N=51; 
cO=0.S$;XN=c0O*xones(1,length(D); 9%dc 
forn=1l:2:N 2% even harmonics are zero 

H=abs(1/(1+j*RC*wOxn)); 
phi=angle(1/(1+j*RC#wOx#n)); 
XN=XN+H#cos(W0O*nxt+phi)*sinc(n#0.5); 
end 
plot(LxN); 
Xlabel([time RC=',num2str(RC)]);grid; 
Set(scaixXtckE35E3 三 机 0 和 瑟 |) 





二 、 周 期 信号 通过 系统 的 响应 





例 2 周期 方 波 通过 RC 系 统 的 响应 。 





-5 -3 
了 
号 0.5 冲 (| 
0 
本 E 





1 
t( 秒 ) 


3 六 
RC=0.1 
1 3 和 

















1 
t( 秒 ) 


三 、 离 散 系统 频率 响应 的 计算 








(ei2 ) = 证 人 
4(ei2) 





计算 频率 响应 的 MATLAB 函 数 
b 分 子 的 系数 a 分 母系 数 


w 抽样 的 频率 点 (至 少 2 点 )， w 在 0~27 之 间 
幅 频 特性 : abs， 相 频 特性 : angle 








三 、 离 散 系统 频率 响应 的 计算 


























例 : 画 厅 (ei2)= 上 的 幅度 响应 曲线 
1--we 1 

b=[]]; 
al=[] -0.9]; a2=[1 0.9]; 
w=linspace(0,2*pi,512); 
h1=freqz(b,al,w); 
h2=freqz(b,a2,w); 
plot(w/pi,abs(h1),w/pi,abs(h2),:); 
legend(\alpha=0.9,\alpha=-0.9); 


三 、 离 散 系统 频率 响应 的 计算 











例 : 画 瓦 (ei2) = 


一 二 的 幅度 响应 曲线 
忆 


1 一 C 





”= 
站 利用 MATLAB 进 行 连续 系统 的 
复 频 域 分 析 


4 部 分 分 式 展 开 的 MATLAB 实 现 
4 9) 的 零 极 点 与 系统 特性 的 MATLAB 计 算 


一 、 部 分 分 式 展开 的 MATLAB 实 现 








Tesidue 


numden 分 别 为 XCs) 分 子 多 项 式 和 分 母 多 项 
式 的 系数 向 量 。 


[为 部 分 分 式 的 系数 ，p 为 极点 ，K 为 多 项 式 
的 系数 。 若 为 真 分 式 ， 则 kk 为 零 。 





二 、 太 (9) 的 零 极点 与 系统 特性 的 MATLAB 计 算 








计算 多 项 式 根 roots 的 函数 可 用 于 计算 尽 (9) 的 零 极点 。 


刀 (9) 零 极点 分 布 图 可 用 pzmap 函 数 画 出 ， 
调用 形式 为 





表示 画 出 sys 所 描述 系统 的 零 极 点 图 。 





例 I 用 部 分 分 式 展开 污 求 KB 的 反 变 换 ， 
将 


2%program7_1 

format rat 多 将 结果 数据 以 分 数 的 形式 输出 
num=[1l 2];: den=[1 430]; 
[:p]=residue(Cnum,den) 


运行 结果 为 =-1/6 ,-1/2 ,2/3 
二 三 3 天 二 由 且 个 
2/3 -0.5 -1/6 
故 F(s) 可 展开 为 X()= 一 一 + 本 
8 S+1 +3 








国 本 
xD= 民 XGOI=3AD-7e VD-Ee MO 


例 2 用 部 分 分 式 展开 法 求 Xs) 的 反 变 换 。 
252 + 3 本 


天 (人 = 一 全 一 一 
E.， (CS+D(02 +S+2) 


和 % program7_2 

num=[2 3 0 5]; 

den=conv([1 1],[L 1 2]); 

% 将 因子 相 乘 的 形式 转换 成 多 项 式 的 形式 
[PK]=residueCnum,den) 

magr=abs(GD) 狗 求 r 的 模 

angr=angle(r) 多 求 r 的 相 角 


例 2 用 部 分 分 式 展 开 法 求 Xe) 的 反 变 换 。 
Xkf = 2 二 王国 
(CS+D(s +S+2) 


T=-2.0000 + 1.1339i, -2.0000 - 1.1339i, 3.0000 
P =-0.5000 + 1.3229i, -0.5000 -- 1.3229i, -1.0000 
ES=2 

Imagr =2.299, 2.2991, 3.0000 

angr =2.6258, -2.6258, 0 


故 Fs) 可 展开 为 
X(8)=2+ 





3 更 2.2991e 26238 2.2991ei26258 
S+1 5S+0.5+jl.3229 ss+0.5-j1.3229 


xD=2600+3e xD 上 +1.1495e "> cos(1.32297 十 2.6258)x(D) 








] 
例 3” 试 画 出 系统 CO0)= Ta 的 零 极 点 分 布 图 ， 
求 其 单位 冲 激 响应 MD 和 频率 响应 态 j 四 ， 并 判断 系统 是 
否 稳 定 。 





num=[1]:den=[1 2 2 1]; 
SyS=tft(num,den); 
poles=roots(den) 
figure(1);pzmap(SyS); 
好 0:0.02:10; 
h=impulsenum,den,b; 
figure(2);plot(bb) 
title(Impulse Respone) 
[Hw]=freqsnum,den); 
figure(3);plot(w,abs(HD) 
xlabel(\omega) 
title(CMagnitude Respone) 


























”= 
站 利用 MATLAB 进 行 离散 系统 的 
Z 域 分 析 


# 部 分 分 式 展开 的 MATLAB 实 现 
s 末 (z) 的 零 极 点 与 系统 特性 的 MATLAB 计 算 


一 、 部 分 分 式 展开 的 MATLAB 实 现 








num,den 分 别 为 XCz) 分 子 多 项 式 和 分 母 多 项 式 
的 系数 向 量 。 

I 为 部 分 分 式 的 系数 ，p 为 极点 ，k 为 多 项 式 
的 系数 。 若 为 真 分 式 ， 则 K 为 零 。 





二 、 克 (z) 的 零 极点 与 系统 特性 的 MATLAB 计 算 








利用 tf2zp 函 数 计算 扩 a) 的 零 极点 ， 调 用 形式 为 
b 和 a 分 别 为 所 az) 分 子 多 项 式 和 分 母 多 项 式 的 系数 向 量 。 
返回 值 为 零点 、p 为 极点 、k 为 增益 常数 。 
互相) 零 极点 分 布 图 可 用 zplane 函 数 画 出 ， 
调用 形式 为 ee 








例 1 将 XG) 用 部 分 分 式 展 开 。 


18 
中 将 < 
多 1 





2%programg8_1 
num=[18]; den=[18 3 -4 -1]; 
[pk] = residueznum,den) 
运行 结果 为 
T =0.3600 ，0.2400 ，0.4000 


=0s000E0.35333 03333 
k = 


故 XCz) 可 展开 为 


0.30 0.24 0.4 
X(Cz) 三 二 史 让 二 区 区 
ESzaa SS3353< LEE SS3353< 





ZE Hz 








侈 全 覃 由 系统 “ ?100520005c+032 俐 全 志 





了 





求 其 单位 冲 激 响应 扣 和 和 频率 响应 Fei9 。 


和 program 8_2 

b =[1 2 1];a =[1 -0.5 -0.005 0.3]; 
figure(]);zplane(b,a); 
num=[0 12 1]; 

den=[1 -0.3 -0.005 0.3]; 
h=impzOum,den); 
figure(2);stem(h) 
Xlabel(Ck) 

title(Impulse Respone ) 
[H,w]=freqzCnum,den); 
figure(3);plot(w/pi,abs(HD)) 
Xlabel(Frequency \omegay) 
title(Magnitude Respone ) 





Impulse Respone 








运行 结果 。 


加 1 Usonoooooooooooaon 





Magnitude Respone 


























中 
0 01 02 03 04 05 06 07 08 09 1 
Frequency 。” 


”= 
业 ，，，，， 了 系统 状态 
量 分 析 


* 微分 方程 到 状态 方程 的 转换 

* 状态 方程 系统 函数 珑 阵 吾 Cs) 的 计算 
4 MATLAB 求 解 连 续 系统 的 状态 方程 
4* MATLAB 求 解 离散 系统 的 状态 方程 


一 、 微 分 方程 到 状态 方程 的 转换 








num,den 分 别 表示 系统 函数 瓦 Cs) 的 分 子 
和 分 母 多 项 式 ; 


A,B,C,.D 分 别 为 状态 方程 的 和 矩阵。 





二 、 状 态 方程 系统 函数 矩阵 瓦 (9) 的 计算 








A,B,C,D 分 别 表示 状态 方程 的 和 阵 。 

开 表示 函数 ss2 引 计算 的 与 第 K 个 输入 相关 的 系统 
函数 ， 即 豆 (9) 的 第 k 列 。 

num 表示 玖 (9) 第 K 列 的 六 个 元 素 的 分 子 多 项 式 

den 表示 再 (9) 公 共 的 分 母 多 项 式 。 





三 、MATLAB 求 解 连续 系统 的 状态 方程 








获得 连续 系统 状态 方程 的 计算 机 表示 模型 
Sys = SS(A,.B,C.D) 
求解 状态 方程 [yto,q]=lsim(sys,X,tq0) 


sys 由 函数 ss 构造 的 状态 方程 模型 

t ， 需 计算 的 输出 样本 点 ,t=0:dt:Tfinal 
x(k) 系统 第 k 个 输入 在 t 上 的 抽样 值 
qd0 系统 的 初始 状态 (可 缺 省 ) 

yCGk) 系统 的 第 k 个 输出 
to 实际 计算 时 所 用 的 样本 点 ; qd 系统 的 状态 








四 、MATLAB 求 解 离 散 系统 的 状态 方程 








获得 离散 状态 方程 的 计算 机 表示 模型 


sys 由 函数 ss 构造 的 状态 方程 模型 

x(:,) 系统 第 k 个 输入 序列 

qd0 系统 的 初始 状态 (可 缺 省 ) 

y(GC) 系统 的 第 k 个 输出 

n 序列 的 下 标 ; qd 系统 的 状态 


或 直接 利用 





例 1 写 出 系统 y (D+57(D)+10y(D=xD 的 





由 [A.B,C,.D]=tf2ss([1],[1 5 10]) 
但 =3 =10 1 
可 得 4 虽 s-|| C=[0 1] D=0 


ER， 


yO=[0 中 | 


引 
0 


例 2 已 知 某 连续 系统 的 状态 方程 和 输出 方程 为 
本 二 | 
92(I) 0 -1|19:(0D) 1 01|x(D) 
0 
GD 1L0 = LL OLD 
其 初始 状态 和 输入 分 别 为 
和 四国 
5 人 诈 IEI 汉人 态 


求 该 系统 的 系统 函数 矩阵 吾 C) 和 输出 。 





计算 系统 函数 敌阵 瓦 (9) 


A=[23:0 三 咯 :B=[Oil LO; 
C=[1 1; 0 -1H;D=[1 0; 10]; 
[numl,den1]=ss2tf(A,B,C,D,1) 
[num2,den2]=ss2tf(A,B,C,D,2) 
运行 结果 

numl=l 0 -=! 


denl =1 -1 一 2 
Tamn 公 三 ON 
OO 
en2 一 二 中 国人 入 
所 以 系统 函数 矩阵 玉 () 为 Po- 加 | 
2 


计算 输出 


%Program 9_1 

A=[2 3;0 -1];B=[0 1; 1 0]; 
C=[1 1; 0 -1];D=[1 0; 1 0]; 
qd0=[2 -1]; 
ddt2: 
X(:,1)=ones(length(tb),1); 
X(::2)=exp(=3#t); 
SyS=SS(A,B,C,D); 
y=lsim(SyS,x,tq0); 
Subplot(2,1,1);plot(by(,1),r); 
ylabel(Cy1l1(0);xlabel(Ct); 
Subplot(2,1,2);plot(ty(:,2)); 
ylabel(Cy2(0 );xlabel(Ct); 























例 3 已 知 某 离散 系统 的 状态 方程 和 输出 方程 为 
ak+I] [0 Iwg] [of 
OKE+HU| -2 32[ 晤 | 11 

人 避 国 刘 
y[K] 祁 本 三川 凡 吕 | 


其 初始 状态 和 输入 分 别 为 


[0] 1 
三 天 三 天 
长 x[K]=x[] 


求 该 系统 的 输出 。 





计算 输出 


%Program 9_2 

A=[0 1; -2 3];B=[0;1]; 
C=[1 1; 2 -1];D=zeros(2,1); 
q0=[1; -]]; 
N=10:;:x=ones(1,N); 
y=dlsim(A,B,C,D,x,q0); 
Subplot(2,1,1);y1=y(:,1) 
Stem((O:N-1),yl); 
Xlabel(k);ylabel(Cyl ); 
Subplot(2,1,2);y2=y(:,2) 3; 
stem((O:N 一 1),y2); 
Xlabel(k):ylabel(y2 7); 
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